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We consider perturbations of a Schwarzschild black hole that can be of both even and odd parity, 
keeping terms up to second order in perturbation theory, for the 1 = 2 axisymmetric case. We 
develop explicit formulae for the evolution equations and radiated energies and waveforms using the 
Regge-Wheeler-Zerilli approach. This formulation is useful, for instance, for the treatment in the 
"close limit approximation" of the collision of counterrotating black holes. 



There has been a recent surge of interest in black hole perturbation theory, stemming from the successful applica- 
bility of such a formalism in the case of black hole collisions (see [|] for a review). In particular, the use of second 
order perturbative computations |^ has proven quite useful to endow the perturbative formalism with "error bars" 
and also to increase its accuracy. 

The computations associated with the second order formulation are quite laborious. In particular, the second order 
perturbative equations are linear differential equations that contain "source" terms that are quadratic in the first order 
perturbations. This forbids us from doing a "generic" second order computation, since it would involve an infinite 
number of source terms. Computations have to be restricted to certain multipolar orders and parities to achieve a 
computable expression. The expressions we have presented in the past j^] are all for even parity axisymmetric i — 2 
first order perturbations. This is a significant case, since for first order perturbations such modes are the ones that 
dominate the initial data for black hole collisions (at least in the head-on case, in the non-head on case one also 
has non-axisymmetric m = ±2 modes). If the black holes individually have spin, however, odd parity terms appear 
as well. If one wishes to study such collisions with the second order formulation one needs to lay out the second 
order equations for odd and even parity perturbations. Since the "source" terms that appear in the second order 
equations are quadratic, one has quadratic contributions from odd parity first order modes in the even parity second 
order equations. The intention of this paper is to present a comprehensive account of the perturbative equations 
involving odd terms, up to second order. In a separate publication we will apply the formalism to the collision of 
counterrotating black holes. 



In order to fix notation and to set the computational philosophy for the rest of the paper, we recollect here several 
results of first order perturbation theory that are known, but that it is useful to have laid out in a compact fashion 
for even and odd perturbations together. 

We consider the axially symmetric perturbations of a spherically symmetric background, up to and including 
second order, using the Regge-Wheeler framework. In first order perturbations there is no loss of generality in 
the consideration of axisymmetric perturbations since modes with different m value decouple. This is not the case, 
however, in second order theory and our treatment is therefore restricted to the axisymmetric case. In terms of 
the usual Schwarzschild-like coordinates t, r, 0, cj), any i = 2 axially symmetric perturbation may be written in the 
"Regge-Wheeler form" |^ 



I. INTRODUCTION 



II. PERTURBATIVE FORMALISM: REVIEW OF SOME FIRST ORDER RESULTS 



A. The Regge-Wheeler decomposition and the Regge-Wheeler gauge 




(1) 



(2) 



(3) 



1 



9te 
9t,p 
gr(i> 
9re 

gee = r 



eVi\r, t) + e^V^^\r, t)\ dP2{0)/d0) 
eh^^Jir, t) + e2/i(2)(^, t)] sm{d)dP2{e) / 86) 
eh['J {r, t) + e^/ifj (r, t)] sm{e)dP2 {6)/ 86) 
eh!^^\r,t)+e''h^^\T,t)] 8P2{6)/86) 



{ 1 + [eX(i) (r, t) + £2^(2) (^r, t)]P2 [6) 
+ [eG(i) (r, t) + £2(5(2) (^^ <)]a2p^ /5^2) I 

|sin2 6l + sin2 6l[e^(i)(r, i) + e^^'^) ^)]p2(6') 
+ sin(6l) 008(6*) [eG^i) (r, i) + e^G^^) (^^ ^)]gp2 (0) /86) } 



560 



e;^^'j(r,t) + t'h\^>{r,t) [cos{d)8P2{e) / 89 - am{6)8'P2i6)/86' 



.21(2) 



(4) 
(5) 
(6) 
(7) 

(8) 

(9) 
(10) 



where P2{6) = 3cos2(0)/2 — 1/2. We introduce the expansion parameter e, which defines the perturbation order, and 
identify the corresponding metric coefhcients with the superscripts (1), (2) and use the subscript "o" to refer to the 
odd parity quantities. 

As we stated above, any axisymmetric perturbation can be put into the Regge-Wheeler form. However, for 
calculational simpHcity, it is usuaUy convenient to take advantage of the coordinate freedom to restrict somewhat the 
form of the metric. An example of such choice is the so called Regge - Wheeler gauge (RWG), where the nonvanishing 
metric coefficients are 



gtt 



Qrr 



-(1 - 2M/r) {l - [eH^^\r,t)+e^H^^\r,t)]P2{e)] 
(1 - 2M/r)-i {l + [eH^^\T,t)+e^H^^\T,t)]P2{e)] 



grt^[eH^^\r.t) + e'H'^>{T,t)]P2{6) 
gte = [e/i[,'j(r,i) +e2/,(2)(^,t)]sin(^)9P2(e)/a0 
gre = [eht^{r,t) + e'hf^ir.t)] sm{e)8P2{9)/8d 



2rr(2)^ 



gee 
94>4> 



{ 1 + [ei^ (r, t) + e^/C (2) ^^jp^ | 

|sin2 6 + sin^ 6[eK'^^\r,t) + e'^K^^\r,t)]P2{6)} 



(11) 
(12) 

(13) 
(14) 
(15) 
(16) 

(17) 



that is, we are demanding that /iq ' — h^^^ = G'-'-' = /ijo = 0, i ~ 1, 2. 

A remarkable aspect of the Regge-Wheeler gauge (RWG) is that it is unique, in the sense that the metric perturba- 
tion coefficients in the RWG can be uniquely recovered from those in an arbitrary gauge. That is, given an arbitrary 
metric, in an arbitrary gauge, there is a well defined procedure that brings it to the RWG. 

Let us exhibit this property for first order gauge transformations. One can uniquely choose a gauge transformation 
vector that brings the metric to the RWG form |e|. The resulting transformation formulae are. 
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Something that should be stressed is that given the uniqueness of the RWG, quantities computed in such a gauge 
(hke, for instance, metric components), are gauge invariant. That is, substituting the above formulae in the definition 
of any quantity in the RWG one obtains expressions for the quantity in any gauge. 

The procedure can in principle be repeated at every order in perturbation theory. To achieve this, one first 
performs a first order transformation that yields the metric (at first order) in the RWG. This transformation will 
induce transformations at all higher orders as well. One then makes a purely second order gauge transformation that 
brings the second order portion of the metric to the RWG. This will leave the first order piece in RWG, and will 
modify the third and higher orders as well. One can continue this procedure up to an arbitrarily high order and the 
metric will be, up to that order, in the RWG. Unfortunately, the procedure is not quite unique at higher orders. The 
reason for this is that if before performing the second order gauge transformation, one performs an ^ = first order 
gauge transformation, the procedure we outlined yields at the end of the day a different second order metric, but still 
in the RWG. That is, the metric one obtains at the end of the procedure is not unique, in spite of being in the RWG. 
Fortunately, one can isolate second order quantities that are invariant under these ^ = transformations. This point 
has been discussed in reference Q in detail, so we will not repeat it here. From the expressions for waveforms given 
in this paper, one can apply the techniques of reference in order to compare with numerical results without the 
^ = ambiguity. 

From now on, we will do computations in the RWG. Because of what we just discussed, there is no issue of gauge 
invariance in our formalism, since all quantities in the RWG (appropriately modified to take into account the £ = 
issue) will be "gauge invariant" in the sense described above. 



B. The first order Zerilli equation 



As shown by Zerilli, we may introduce a function ^p'^^^{r,t), such that if we write for the metric components in the 
RWG, 
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then the linearized Einstein equations are satisfied, provided only that if)'^^^ (r, t) is a solution of the Zerilli equation 
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It is straightforward to obtain "inversion" formulas for ?/''^^'(r, t) in terms of the metric coefficients. From ( p^ and 
(p6|), we have 
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while from (E4) and (E9) we find 
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3(2r + 3M) V dr 

Equations pT| ) and (^) are equivalent, as far as first order perturbations are concerned, up to an additive function 
of r. This ambiguity is of no concern since the relevant physical quantities, like radiated energies and waveforms are 
all computed in terms of the time derivative of ?/'^^\ so it is conventionally ignored, we will assume we set the additive 
function to zero. Equation ( ^2[ ) can be used to obtain a gauge invariant form for ip^^^r^t), i.e., an expression for 
ip''^'^{r,t), valid in an arbitrary gauge. This is an especially important result, because it allows us to obtain the initial 
data for the Zerilli equation, directly in a n ap propriate gauge. It should be noted, however, that although ip^^') (r, t) 
is uniquely defined in a general gauge by (^2|), we do not have, in a general gauge, unique expressions for the metric 
perturbation functions in terms of 'il)'^^\r,t). The general statement of Zerilli's results is, in this case, that if the 
metric perturbations satisfy the first order Einstein equations, then ip^^^r^t) satisfies the Zerilli equation (|2^). 

The situation, as we shall see, is somewhat different for the second order perturbations, because of the presence of 
"source terms" in the Zerilli equation. 



C. The Odd parity equation 

One can also write a "Zerilli" -like equation (known as the Regge- Wheeler equation, historically it was discovered 
earlier than the even parity one ^ ) for the odd part of first order metric perturbations and the equations they satisfy. 
If we define 
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the Einstein equations are satisfied, provided that (5'^^^(r, t) is a solution of the equation 
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It is easier to obtain inversion formulas in this case that in the even parity case. We can write, 

Q...,,,. 4 (3a, 
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We can make the same comments as in the even case. Equation (|3^ ) can be used to give a gauge invariant form for 
the odd parity Zerilli function. Both equations allow us to give initial conditions for Q''^^{r, t) in an any gauge. As in 
the even case, there is no unique way to write the metric perturbations in terms of Q^^\r^ t) in an arbitrary gauge. 



III. SECOND ORDER PERTURBATIONS 

We define the second order Regge-Wheeler gauge by imposing that — h\ — G^''' — /120 = 0, i = 1, 2. It is easy 
to check that this gauge is uniquely defined, as long as we restrict to i = 2 axisymmetric first and second order gauge 
transformations, but the general expressions for the second order gauge transformations are rather lengthy. 
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It should be clear from the perturbation expansion, that the second order Einstein equations imply, in turn, that 
the second order perturbation functions satisfy the same linear set of equations as the first order perturbations, but 
with additional terms quadratic in the first order perturbations, that may be thought of as representing "sources" for 
these equations. 

Therefore if motivated by the corresponding expression for the time derivatives dip''^'' /dt, dQ^^'^ /dt, of the Zerilli 
functions in first order perturbation theory, we introduce gauge invariant functions, given by 



2r + 3M 



dt 



r J 



t) 



(40) 
(41) 



it follows that if the second order Einstein equations are satisfied, then ^(^^(r, t), 0*^^^ satisfy equations of the form 
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where Sz, Sq are quadratic polynomials in the first order functions and their t and r derivatives, with r dependent 
coefficients. Sz,Sq may be considered as a kind of "source terms" for the otherwise homogeneous equations satisfied 
byx(2)^0(2). 

But, precisely because of the presence of these "source terms", if we redefine x(2)(r, i), 9(2) (r,i) by the addition of 
a quadratic polynomial in the first order functions and their derivatives, the new, redefined x'^H'^iO; ©(2^'":^) will 
still satisfy an equation of the form (^), provided only that the first and second order Einstein equations are satisfied. 

This arbitrariness in the definition of the Zerilli functions associated to the second order perturbations may be used 
to simplify the form, and asymptotic behavior of the "source terms" Sz, Sq. As was discussed in [| , the source terms 
are delicate in the sense that they can be divergent for large values of t and r. The divergence of the source does not 
imply a physical divergence in the problem, but it poses serious difficulties for numeric computations. If we make the 
following choice, based on the study of the solutions of the Zerilli equations for large values of r, the source simplifies 
and also is well behaved at spatial infinity. 
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as our definition for x(2''(?', t), while 9(2) (r,i) is defined by 
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where ■0 = '0''"'^'', Q = Q'^-* and the subindices indicate partial derivatives. 

IV. GRAVITATIONAL WAVE AMPLITUDE 

A straightforward computational way to define a gravitational wave amplitude, is to consider a gauge in which the 
space-time is manifestly asymptotically flat. For instance, consider the typical gauge discussed in reference ||^. In 
this gauge the metric has components that behave in the following way, 

H2 - 0(l/r3) 

h, ~ 0(l/r) 

K ~ 0(l/r) 

G~0{l/r) 
h,o ^ 0(l/r) 

h2o ^ Oir) (47) 

and, of course, Hq = Hi ^ ho — h^o — 0. The gravitational wave amplitude can then simply be read off from the 
components K , G and /i2o- We therefore need to relate the behavior of these components of the metric to the Zerilli 
function in an asymptotic region far away from the source. 
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In order to study the asymptotic behavior, one considers a solution of the Zerilh equation in terms of an expansion 
in inverse powers of r with coefficients that are functions of the retarded time t — r* , as discussed in detail in 
Substituting in the Zerilli equations we obtain, 
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For this behavior of first order Zerilli functions, the sources Sz, Sq for second order Zerilli functions x'^^ 0'^' 
decay as 0(l/r^) asymptotically. Therefore, x^^^G^^-* will behave in a similar way to V'^"'^'') Q^^^ for large r: 
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We can now compute the asymptotic form for the metric perturbations, to first order and second order, in Regge- 
Wheeler gauge, using ( |4^ ) and (^) in ( p5| ) and (^^. Unfortunately, in the Regge- Wheeler gauge the metric is not 
in the manifestly asymptotically flat form (47). We therefore need to perform a gauge transformation to bring it to 
such a form, to first and second order. 

In order to do this, we start by recalling that for any £ — 2 axisymmetric gauge transformation, we can write the 
gauge transformation vectors as. 
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To achieve the asymptotically flat gauge ( ^ ) we demand that H^^\r,t) = H['^\r,t) = h^o\r,t) = h^^ 
the general form of the gauge transformation equations we have, 
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where the quantities with tildes are given in the Regge- Wheeler gauge, {M^^\ M^q\ M^i"^ } , are the components of 
the even parity gauge transformation vector, and Ad'^^ that of the odd parity gauge transformation vector. 

We can see that, given iJg^^ and h[^\ the first and second equation above may be used to solve for M^^^ and 

(1) 



A4\ ' . Then, the third equation leads to a partial differential equation for A4^^^ while the equation for h^J fixes M. 
Therefore, the set of equations, and therefore the gauge choice, is consistent. The solution, assuming only that Hq^' , 
h[^\ h'^\ and ^(1) are given, is non unique since one has to integrate partial differential equations, so the solution 
in general involves free functions. One may choose such functions in a judicious way to ensure that the non-vanishing 
metric components behave in a manner consistent with asymptotic flatness. 

Similarly, wo define the second order asymptotically flat gauge by imposing the flrst order asymptotically flat gauge 
conditions and H^ '{r,t) = 0, Hf'{r,t) = 0, h}y '{r,t) = and h};J{r,t) = 0. Clearly, the same reasoning as above, 
applied to the second order gauge transformations, shows that this is a consistent choice. 

The end result of these calculations is that one may write the first order gauge transformation vectors in the form, 
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The second order gauge equations are more complex then the first order ones. It can be seen that the equations can 
also be solved itcrativcly for all orders in r. but for simplicity we only compute the terms that are going to be relevant 
in the gravitational waveforms. The meaningful terms arc given by 
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and if we compute them in terms of F^,Fq, we finally get 
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We can now compute the gravitational waveforms by reading off the appropriate metric components. For first order 
we find, 
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(iv). 



dt 

-2r{Fl\t-r*) + 



15^ 



and in terms of first and second order ZeriUi functions, we get 



dt 



dt 



0(l/r" 



dK(^\r,t) _ 9G(i)(r,t) 
— o 



dt 

dh^2Jir,t) 
dt 



-2r 



dt 

dQ^^Hr, t) 
dt 



+ 0(r") 



and 



dt 

dK(^\r, t) 
dt 

dh'i}{r,t) 
dt 



^ 7dt 



5^(1) 
dt 



at 



-+0(l/r2) 
r 



^ 3 



dt 



+ 0(l/r2) 



(56) 



(57) 



(58) 



(59) 



It is straightforward to compute the radiated energy, since we are in an expHcitly asymptoticaUy flat gauge. One 
can therefore simply apply the formulae stemming from Landau and Lifshitz pj7|, 



dPower 



1 



1/ 

dt 



dfl 167rr2 
For the perturbations we are considering, 

hgg = {[eX(i)(r,<) + e^K'^^\r, t)]P2{e) 

+ [eG(i)(r, t)+e^&^\r, t)]d^ P2i9) / d9^)^ 



^ 1 / a , 1 5 , ^ ^ 

-i \dt sm- 



^^edf 



(60) 



(61) 
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i^^ = {sin2 0[ei^(i)(r, t) + e^K'^^\r, t)]P2{e) 

+ sin(6') cos(6l) [eG(i) (r, t) + e^Gls) (^^ i)]aP2 (6')/a6i) } 



1 r 



,21(2), 



(62) 
(63) 



and substituting, we get, 
Power 



+ 



3 

10 

36 
35 



' at 



e(2) 



ld_ 

7di 

1 

" 7 



9i ; 7 V 5i 




9i at 



(64) 



As can be seen, the energy has two quadratic contributions, one coming from the hg^ component of the metric and one 
from the diagonal elements of it. They correspond to the "x" and "+" polarization modes of the gravitational field 
and therefore we can consider the expressions in the square brackets as the "waveforms" of the gravitational waves 
emitted. Contrary to the first order case, in which the waveform is directly given by the Zerilli (or Regge- Wheeler) 
function, for second order corrections the relationship is more involved. 



V. INITIAL DATA 



Initial data for solutions of Einstein equations is usually given in terms of the initial value of the three-metric and 
the initial value of the extrinsic curvature of the initial Cauchy surface. For the case of initial data for colliding black 
holes, the popular families of initial data of Bowen and York j|] and punctures |^ naturally come cast in a gauge in 
which|^ Hq = 0, Hi = 0, ho = Q and /iqo = 0. We will therefore present the appropriate formulae that would allow us, 
given a metric and extrinsic curvature cast in such a gauge, to provide the initial data for the first and second order 
Zerilli functions and their time derivatives. 

The extrinsic curvature Kab of a spatial slice is defined by 



Kab = n(c;<i) ('^'^Q + n^Ua) {6"^ b + n'^m) 



(65) 



where is the unit normal to the Cauchy surface where the data are specified. One can relate the extrinsic curvature 
to the time derivative of the three-metric through, 



dt 



= -2NK, 



(66) 



That is, if through some procedure we are given the extrinsic curvature and the three metric, we can compute the 
time derivative of the three metric. Re-expressing this in terms of components of the tensor spherical harmonic 
decomposition due to Regge and Wheeler we get. 



dh']^^/dt ^ 
dh['J/dt = 



-2^/1 - 2M/rK[l\^ dH^'^/dt = -2^1 - 2M/rK^ 
-2^1 - 2M/rK[^^ dG^^ydt ^ -2^1 - 2M/rK^ 



-2^1-2M/rKlll.^ dh^^'J/dt = -2^1-2M/rK^ 



(1) 



(1) 

(h2o) 



(67) 



for first order components. are spherical tensor harmonic components of the first order extrinsic curvature 

tensor, written in the Regge-Wheeler notation that we used in formulae (p}]lO|). For second order we get. 



dhf^/dt = 



-2^1 - 2M/rKf^^^ dH^^^/dt = 



r(2) , 



-2^1-2M/rK\^^ dG^^^dt = 
-2^1-2M/rK^^l^ dh^i^/dt = 



-2^ 1-2M/r Kj^^^ 
-2^1 - 2M/rK^" 



-2^1 - 2M/rK^ 



(2) 
(2) 

{h2a) 



(68) 



^ This happens to coincide exactly with the same gauge condition of the asymptotically flat gauge. 
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(2) 

Here K^J^^ are spherical harmonic components of the second order extrinsic curvature tensor. Note that the expressions 
are similar to the first order ones only in the gauge we are considering. In other gauges (in which the perturbations 
of lapse and shift may not vanish), there will be extra terms, quadratic in first order elements. 

A case in which the formalism just developed will be of use (and we will discuss in the forthcoming publication) 
is that of the collision of counterrotating Bowen-York holes (the cosmic screw case). In such a case, to first order 
the metric has even portions only (similar to those in the Misner initial data) and the extrinsic curvature odd 
portions. So putting together the above formulae, for this case one gets, 



/(2) 



dt 



9(2) 

dt 



2 y/r - 2M 
^7(2r + 3Af)V^ 

3 (if(i))'Af 



-r + 2 M 



7rK, 



(2) 
(hi) 



2 U« 



-r + 2M 



UkI^.M + 21 rK^^^M 



1 r-2M 
14r4(2r + 3M) 



140 A'/r/if ^ + 84 A/V^G*^) + (kI^A - 15 (k^^A - 2 A/r* (k^^ 



2 Mr^ ( K^^ 



96 r if, 



'(1) 

"(h2o 



276 M K, 



'(1) 



144 Af^r K, 



41) 



[hla) (h2a) 



28 r'Gil^ - 28 r4G(2) - UAMrK[ll^K[ll^ - 24 Mrif«^)X«^)^+ 

84M^r^Gi^J + UMr^Gf^ - 42 Afr^G^.^' - 168 A/V/if,^, - 28 Mr^hfl + 56 r^/if^- 



168 M^h^^^ - 56 r^h\"' - 28 r-'H^"' - 28 Mr^H. 
i-r + 2 Af ) (-2 h^iJ + rh'il + 2 r/^fj) 



.2^(2) 



3ff(2) 



r(2) 



7 r- 



2M 



2r3 



^2t^{1) K^i) ^-iAr^frC^) _ 7^2^(2) 



35 K 



(2) 

(/l2o) 



M + 14 rK, 



(2) 



(^^^r'^^lt)^! + 8rK^'^K[ll^^M + UrK^^l^^M + 28 rifgJ^^M ~ 4 r^if^^X^^)^ 
3r2if(i)xW , +8rA'«if,(,^^ 20if«x(i) ,Af-2rifW ^(i)^,^ 

' ('l2o) ('i2o) {h2o) (hlo) 



(69) 



VI. CONCLUSIONS 



We have written explicitly the second order perturbative equations for axisymmetric £ — 2 perturbations including 
odd parity modes. We set up explicitly the evolution equations and the formulae for generating initial data starting 
from a metric and an extrinsic curvature. We also performed an asymptotic analysis to obtain formulas for radiated 
energies and waveforms. These formulae will be useful for considering perturbatively the evolution of spacetimes 
arising from the collision of nearby black holes. They might also be useful for other comparison with nonlinear 
evolutions, as the ones considered in mU. 
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